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In the Horndeski's most general scalar-tensor theories, we derive the three-point correlation 
. . . function of scalar non-Gaussianities generated during single-field inflation in the presence of slow- 

ed ' variation corrections to the leading-order term. Unlike previous works, the resulting bispectrum is 

valid for any shape of non-Gaussianities. In the squeezed limit, for example, this gives rise to the 
, same consistency relation as that derived by Maldacena in standard single-field slow-roll inflation. 

04 ' We estimate the shape close to the squeezed one at which the effect of the term inversely proportional 

^ ' to the scalar propagation speed squared begins to contribute to the bispectrum. We also show that 

^ , the leading-order bispectrum can be expressed by the linear combination of two convenient bases 

whose shapes are highly correlated with equilateral and orthogonal types respectively. We present 
concrete models in which the orthogonal and enfolded shapes can dominate over the equilateral one. 
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I. INTRODUCTION 



The potential presence of primordial non-Gaussianities in the CMB temperature anisotropics can be a powerful 
Oh' probe for the physics in the early Universe-especially for inflation [l[. The inflationary paradigm generally predicts 
^ ■ nearly scale-invariant density perturbations with a suppressed tensor-to-scalar ratio, whose prediction is consistent 
I . with the CMB power spectrum measured by COBE Q and WMAP The detection of scalar non-Gaussianities 

not only breaks the degeneracy among many inflationary models, but it also offers the possibility to discriminate 
O^l , between the inflationary paradigm and other alternative scenarios (such as curvaton Q) 

^ ' There are several different shapes of non-Gaussianities depending on the wave numbers fci, A;2, and satisfying 
the condition fci -I- fc2 + ^3 =0 [s^-jl^. The simplest one is the so-called local shape, which has a peak in the 
squeezed limit (i.e., the limit where the modulus of the momenta approaches fca — > and fci ~ ^2). The second shape 
' corresponds to the equilateral configuration with a peak at fci = ^2 = ^3- A factorizable shape whose scalar product 
with the equilateral template vanishes is called the orthogonal one. There is another shape dubbed the enfolded one, 
which is a linear combination of the equilateral and orthogonal templates. 
, From the bispectrum A-r of the three-point correlation function of curvature perturbations 7?., the non-linear 
parameter characterizing the strength of non-Gaussianities is defined by /nl — (10/3) ^^^j^ fcf . For purely 
adiabatic Gaussian perturbations we have that /nl = 0, but the presence of non-Gaussian perturbations leads to the 
deviation from /nl = 0. The WMAP 9 year data provide the following bounds on the non-linear parameters of local, 
^ \ equilateral, and orthogonal non-Gaussianities, respectively [46j : 
< 

_<^: /ilj°L'^' = 37.2 ± 19.9 (68%CL), /^"l"''' = 37 ± 40 (95%CL), (1) 

/^r'^ 51 ± 136 (68%CL), /^'["" = 51 ± 272 (95%CL), (2) 

/nl""" -245 ± 100 (68%CL), /°'^^° = -245 ± 200 (95%CL). (3) 

Since the non- linear parameter of the enfolded shape is given by /^if"'^ ~ (/nl"' " /nl''°)/2 13 113 j '^e obtain the 
following bounds from Eqs. ^ and ([3]): 

/™f°id ^ -^43 ^ ^gg o;^ , /™f°id = 148 ± 236 (95 % CL) . (4) 

For the local, orthogonal, and enfolded shapes the model with purely Gaussian perturbations (/nl = 0) is outside the 
68 % observational contour, but, apart from the orthogonal case, it is still consistent with the WMAP constraints at 
95 % CL. 

In standard single-field inflation based on a canonical scalar field, Maldacena [l^l showed that the non-linear 
parameter in the squeezed limit is given by f^{f^ = (5/12)(l — tt-k), where n-ji is the scalar spectral index. Creminelli 



and Zaldarriaga [47l | pointed out that the same non-Gaussianity consistency relation holds for any single-field model 
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under the condition that only one mode of curvature perturbations survives after the Hubble radius crossing while 
the other one decays^ (see Refs. jsoMssj for related works). 

In the context of single-field k- inflation (53 |. the bispectrum of curvature perturbations was first derived by Seery 
and Lidsey in 2005 [ll] . Since the scalar propagation speed squared can be much smaller than 1 , it is possible 
to realize the large equilateral non-linear par ameter I/^l"''! ^ V^s S> 1. If we naively take the squeezed limit for the 
leading-order bispectrum derived in Refs. |l8l.[56j. the term proportional to does not disappear. This comes from 
the fact that the slow-variation corrections to the bispectrum need to be taken into account to estimate the local-type 
non-Gaussianity correctly. In fact, Chen et al. [s^ showed that the Maldacena's consistency relation is recovered in 
the squeezed limit by carefully computing all the possible slow-variation corrections to the leading-order bispectrum. 
Thus the slow-variation single-field k-inflation models with ^ 1 lead to small local non-Gaussianities, even though 
the equilateral non-linear parameter can be large. 

In the Horndeski's most general scalar-tensor theories with second-order equations of motion [57l - l60j . the leading- 
order three-point correlation function of curvature perturbations was derived on the quasi de Sitter background [6ll l62j 
(see Refs. [63l - l65| for the scalar non-Gaussianities in related Galileon models and Refs. [g^ for the bispectrum of tensor 
perturbations in the Horndeski's theories). Although the result is valid for the estimation of the equilateral non-linear 
parameter, the bispectrum is not general enough to be used for any shape of non-Gaussianities. In this paper we take 
into account all the possible slow- variation corrections to the leading-order bispectrum in the Horndeski's theories^. 
Not only we reproduce the Maldacena's consistency relation in the squeezed limit, but we identify the shape close to 
the squeezed one at which the term l/cj. begins to contribute to the bispectrum. 

Given our general expression of the bispectrum in the most general single-field scalar-tensor theories, we can 
evaluate the non-linear parameters of several different shapes to confront each inflationary model with observations. 
In particular the result I/n^^^'I ^ 1 is robust for any slow- variation single- field model, so the detection of non- 
Gaussianities in the squeezed limit will allow us to falsify the slow- variation single-field scenario. Note that in realistic 
observations the shape is not completely squeezed, in which case the bispectrum can be affected by the appearance 
of the term l/c^ mentioned above. Our results are useful to distinguish such difference accurately. 

If c2 < 1, then the non-linear parameters 1/^1" |, |/nl^°|, and can be much larger than the order of 

1. Which shape dominates over the other ones depends on the models of inflation. In Ref. [631 '^^s shown that 
the correlation between the equilateral template and the shapes arising from the Horndeski's theories is quite high, 
but linear combinations of equilateral operators can give rise to a significantly different shape for a wide range of 
coefficients [13, . In this regard we anticipate that there may be some models in which the shape orthogonal to the 
equilateral template provides an important contribution to the bispectrum. 

In this paper we show that the leading-order three-point correlation function in the Horndeski's theories can be 
expressed by a linear combination of two bases whose shapes are highly correlated with equilateral and orthogonal 
shapes respectively. This decomposition is useful because the contributions from the equilateral and orthogonal shapes 
can be easily estimated for concrete models of infiation. We show that in k-inflation with the covariant Galileon terms 
there are cases in which the correlations with the orthogonal and enfolded templates are larger than that with the 
equilateral one. Thus the shapes of non-Gaussianities allow us to discriminate such models from observations. 

This paper is organized as follows. In Sec. |ll] we review the background and linear perturbation equations in 
the Horndeski's theories. In Sec. IIIII we derive the three-point correlation function of curvature perturbations in 
the presence of slow- variation corrections to the leading-order bispectrum. In Sec. IIVI the non-linear parameter /nl 
is evaluated in the squeezed, equilateral, and enfolded limits, respectively. In Sec. |V| we express the leading-order 
bispectrum in terms of equilateral and orthogonal bases. In Sec. IVII we show concrete models of inflation in which 
the orthogonal and enfolded shapes can dominate over the equilateral one. Sec. IVIII is devoted to conclusions. In 
Appendix we show the details of the slow-variation corrections to the bispectrum. 

II. EQUATIONS OF MOTION FOR THE BACKGROUND AND LINEAR PERTURBATIONS 

The action corresponding to the most general scalar-tensor theories is given by [13, Ull 



S ^ d XsJ—g 



-^R + P{c^, X) - Gsic^, X) n^ + Ci + £5 



(5) 



^ If the decaying mode is non-negligible relative to the growing mode, the Maldacena's consistency relation can be violated [43. l49l| . 

^ In the effective field theory of inflation (which allows the equations of motion higher than second order) , a similar approach was taken by 
Cheung et al. [50l | to show that the Maldacena's consistency relation holds in the squeezed limit. While the authors in this paper mainly 
focused on the local shape, we derive the full bispectrum in the Horndeski's theory which can be used for any shape of non-Gaussianities. 
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where g is the determinant of the metric g^j^, Afpi is the reduced Planck mass, i? is a Ricci scalar, and 

£4 = G4(0,X)i? + G4.x[(n0)'-(V^V.0)(V^V'' (6) 
£5 - G5(0,X)G^,(V'^V>)-iG5,x[(n<^)'-3(n0)(V^V,0)(V^V>) + 2(V^V«0)(V"V^0)(V'^V^0)]. ^ 

Here P and Gi's (i = 3,4,5) are functions in terms of and X ~ — 9^(/)9^0/2 with the partial derivatives Gi,x = 
dGi/dX, and G^ii/ = R/^^i, — is the Einstein tensor (i?^^ is the Ricci tensor). 

We consider the following ADM metric [G^I with scalar metric perturbations a, ip, and TZ about the flat Friedmann- 
Lemaitre-Robertson- Walker (FLRW) background 

ds^ = -[(1 + af - a{t)-^ e-"^^ [dil^f] df + 2d,il^dtdx' + a{tfe^'^dx^ , (8) 

where a{t) is the scale factor with cosmic time t. We choose the uniform field gauge Scf) = 0, which fixes the time- 
component of a gauge-transformation vector . The spatial part of is fixed by gauging away a perturbation E 
that appears as a form E^ij in the metric (|5]). 

The background equations of motion are given by 

ml^H'^F + P + 6HGi^^<j) + (G3.0 - 12 H^Gi^x + 9 H^G^,^ - Px) 

+ (6 G4.^x - 3 G3,x - 5 Gs.xi?') + 3 (Gg.^x - 2 G^^xx) - ff'Gg.xx^' = , (9) 
(1 - AScix - 2(5g5X + 2(5G50)e = Spx + SSqsx - 2(5g30 + 6 5g4x - Sg44, - 6 Sg54> + 3 Sg5X + 12 Sqaxx + 2 Sqsxx 

— 10 (5g40X + 2 (5g400 ~ 8 <5g50X + 2 (5g500 ~ '5</.('5g3X + 4 (5g4X — '5g40 

-1-8 Sgaxx + 3 5g5X - 4 (5g50 + 2 Sqsxx - 25ga4,x - 4 5Gb<t>x) , (10) 
where H = a/a is the Hubble parameter (a dot represents a derivative with respect to t), F = 1 -I- 2Gi/M^^^ and 

^I'-JT^^ ^'-''^M^' = 
Gi^^cf) GA^cj,x4>X Ga,(I,(I,X Gi^xxX"^ ^ Gs.^X 

'^^^^ = iWp?^^ ' '5g4,X - , ^G,,, = J^T^. Sg4XX= , -^GS^-^^, 

r G^^xHcjiX Gz_xxH4>X'^ G^^^xX'^ c Gs,.^,p4>X 

6g5X = , <5g5XX - ^2 ^ <5g50A- - ^^^p^ , <5g500 = ^^2^^ , (H) 

whose magnitudes are much smaller than 1 during inflation. The terms Scirpx, Scirfxf,, Scbrpx, SG^^^ as well as 
Sg3,PX = Gs^^xX'^/iM^.iH^F) and (5g300 = G3^^^<i}X/{M^^iH^F) are second-order of e. From Eq. ^ it follows that 

e = <^px + 3(5g3X - 2(5g3</. + 6 Sgax ~ <5g4</, - 6 (5g50 + 3 5g5X + 12 Sgaxx + 2 5g5XX + C>{e'^) . (12) 
For the quantity Sp = F/{HF) we have 

Sf = 2(5g40 + O(e') . (13) 

Using the relations between TZ, tjj, and a that follows from Hamiltonian and momentum constraints, the second-order 
action for perturbations reduces to (60l - [6^ Igs} 



/ 



52 = / dtd^xa^^ 



7^2-4 (971)= 



(14) 



where 



wi(4wiW3 -H9w^) 
^ = ' 



2 Z{2w\w2H — W2WA + '^Wi'WiW2 — 2w\w2) 

~ ?z;i(4w;iu;3 9w|) ' ^ 
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and 



wi = M^iF ~4XGi^x ~2HX^G5,x +^XG5,^, (17) 

W2 = 2M^iHF -2X^G3^x'16H{XGi,x + X^Gi^xx) + 2^iG^^^ + 2XGi,4.x) 

- 2H^{5XG5.x + 2X'^G5.xx) + AHX{2,G^.^ + 2XG5.,^x) , (18) 

W3 = ^9M^,iH^F + 3iXPx + 2X^Pxx) + l8H^{2XG3,x + X^G^^xx) - 6X(G3,0 + XG^.^x) 
+18H^{7XGi^x + IGX^Gi^xx + ^X^G^^xxx) - l8H^{Gi^^ + hXG^.^x + 2X^Ga.^xx) 
+ 6H^4>{l^XG5,x + ISX^Gs _xx + 2X^G^zxxx) — 18i7^X(6G5^0 + QXGs^^x + SX^Gs^^xx) , (19) 

Wi = Ml^F -2XG^.^-2XG^^xii}- (20) 

For later convenience we introduce the following parameter 

_ Qcl 



Ml,F 



e + 5g3X + (5040 + 85gaxx + (^GSJf + 25g5XX + C(e ) • (21) 



At linear level the curvature perturbation obeys the equation of motion 



5C2 

We decompose TZ into the Fourier components, as 



at 



TZ(t, x) = ^-i— / d^kTZ(T, fc)e*'=-^ , 7^(T, k) = u(t, k)a(k) + u*(t, ^k)aU-k) 
{2Try J 



(22) 



(23) 



where r = J a ^ dt, k is the comoving wave number, a{k) and a^(fe) are the annihilation and creation opera- 
tors, respectively, satisfying the commutation relations [a{ki) , {k2)] = {2tt)^ S'-^^ki — ^2) and [a(A;i), 0(^2)] = 
[at(fei),at(fc2)] =0. 
Introducing a rescaled field v = zu with z = a^2Q, it follows that 

c^fc^-y)^' = 0, (24) 

where a prime represents a derivative with respect to r. Under the slow-variation approximation the term z" / z can 
be expressed as 

^ = 2(ai7)2 (1 - r + i''-'^ \') + ' ^2^) 

where 

^^^=i?(^^''^ + ^^' '^■^ = ^' '^=^- ('^^ 

Taking the dominant contribution in Eq. (j25p and using the approximate relation a ~ —1/{Ht), we have z"/z ~ 
The solution to Eq. ((24|). which recovers the Bunch-Davies vacuum state (w = e~*^='^'"/-\/2c7fc) in the asymptotic past 
(fcr — > —00), is given by 

^("'')- 2(c.fc)3/^VQ ^^ + "-'"^- ^''^ 

The slow- variation terms in Eq. (P5|) provide the corrections to the mode function (P7|) . Later we shall discuss the 
effect of such corrections on the primordial non-Gaussianities. 

The power spectrum Pn{ki) of curvature perturbations, some time after the Hubble radius crossing, is defined by 
(0|7^(0,fcl)7^(0,fc2)|0) = i2TT'^/kf)Vn{ki) (27r)3j(3)(fei + k2). From Eq. ^ it follows that 

= o 2S F ' (28) 
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which should be evaluated at Cgk = aH . The spectral index n-ji is given by 

d\nVn 



d\nk 



(29) 



The corrections to the solution (f27| only give rise to the O(e^) terms in Eq. (f29| . 

Similarly the power spectrum Vh and the spectral index nt of gravitational waves are given, respectively, by [60l . |62 



if2 



d\nVh 



2n^Qtc: 



d\u.k 



ctk—aH 



where Qt = «;i/4 = (A/2^/4)(l 
both Viz and 7^/i remain constant, the tensor-to-scalar ratio can be evaluated as 

r = — 
Vn 



4(5g4X - 2(5g5X + 2Jg5</.) and = ?«4/wi 



16c,,e,, 



-2e - 5f , 
1 + 45g4x 



(30) 

2(5g5X - 4(5g50- When 
(31) 



III. THREE-POINT CORRELATION FUNCTIONS IN THE PRESENCE OF CORRECTION TERMS 

In the Horndeski's theories the third-order action of perturbations was derived in Ref. [6l|, [G^I- Here we do not 
repeat the details, but we summarize the main results. Under the approximation that all of the slow-variation terms 
in Eq. (|lip are much smaller than 1, the third-order action reads 

^3 = j dtd^xiya^CiMl{Rn^ ^ aC2Ml{R.{dnf + a^C^iM^iil^ 

+ aCeTZ^d^TZ + C7 [d^n{dnf - TZd,dj{diTZ){djTZ)] /a + a{Cs/Mpi) [d'^ndMX - TZdidj{d^TZ){djX)] 

where d^X = QTZ. The dimensionless coefficients C; (i = 1, • • • , 8) and the coefficient J^i are [i^l 

Ci = f - (--1) +^{es-Vs- ^Sg3X - 12<5g4x - 32 + 12(5g50 - lOW - 86g5Xx) + ^(e^), 



(33) 

C2 = FeJ^-l] +^{es + Vs-2s + 45g4X + 2Sg5X - 4(5g50) + 0{e^), (34) 



I \ Fe 



C3 = —5^-77- ( ~ ^ 1 ^ ~ ) + ~jr~S~\ ~\^G3X + 4(5g4X + i&GhX - ^G4<i> " 4(Sg5</> + 8(5G4JfX + ^Eg'oXx) 



Fe,Mpi/^ 2A\ j^e.Mpi f 1 

- (3 + 2A3x) <5g3X - 8 (5 + 2A4x) fc4xx - 4(4 + \'hx)^G^xx 

+<5g40 + 8,5g50 - 8,5g4X - 9W - 6^ [(I + \3x)^l^x + ^('5')] | + O(e^), (35) 

C4--^+0(6^), (36) 
1 



<^5 = 7T;(e. -4<5G3X-8fc4XX + 8fa + 4(5G5Xx) + 0(e'), (37) 
Ce ^ 2f(^] {{l + \3x)SG3X+H3 + 2\4x)SG4XX+SG5X + {5 + 2X5x)SG5Xx]+0{e^), (38) 



■3X — OUG4XX -I- OUG5X T ^ug5XX) T i^yC^'' 

2 

^ H \ 

C7 = -^^(^) iSG3X+(iSG4XX + SG5X+SG5Xx) + Oie^), (39) 

Cs = 2^{dG3X+^SG4Xx)+0{e^), (40) 
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where Xsx = XG^^xx/Gs^x, A 



4X 



XG4^xxx/G4,xx, A5X — XG^^xxx IGz.xx, Mi = ^AfiF — 2A{XGi^x 



X^Gi.xx) - 6H<j>{5XG5,x + 2X^G^,xx) + 6X(3G5,0 + 2XGz^^x), ii = 2u;i/u;2, and 
?i;i(4ti;iti;3 + Owa) 



(42) 



A = 



[ZX^P^xx + 2X'^Pxxx ,xxx) — 2(2X^G'3^0x + X^G'i^^xx) 



+QH^{QX^Gi^xx + ^QX^Gi^xxx + 4X'^G'4_xxxx) — 3H(j){3XGi^^x + 12X^G4,0xx + 4X^G4,0xxx) 
+H^{3XG5,x + 27X^G5^xx + 24X^G5.xxx + 'iX^G^^xxxx) 

^6H^{6X^G5^^x + QX^G^^^xx + 2X''G5^^xxx)] ■ (43) 

The exphcit form of the second-order term £,{S^) in Eq. (|35|) is given in Appendix of Ref. [62j. The coefBcient J^i 
involves the terms with the spatial and time derivatives of TZ and X. These provide the corrections to the three-point 
correlation function higher than first order in slow- variation parameters'^. Since we are interested in the bispectrum 
up to first order, we neglect the contribution of the term Ti{SC2/STV)\i in the following discussion. We also evaluated 
other boundary terms and found that they only lead to the contribution higher than the order e. 
The vacuum expectation value of TZ for the three-point operator in the asymptotic future (t — > 0) is 

(7^(fel)7^(fc2)7^(fc3)) / dra (0| [7^(0,fcl)7^(0,fc2)7^(0,fc3),Hint(T)] |0) . 



(44) 
J —00 

The interacting Hamiltonian Hint is related to the third-order Lagrangian £3 as Hint = ^ £3, where ^3 = f dt C3. We 
write the three-point correlation function in the form 



where 



{n{k,)n{k2)n{ks)) = (27r)35(3)(fci + ^3 + fc3)(7'K)'-^7j(fcl, ^2, fca) , 



^n{ki,k2,k3) = ^3^\ „ -4R(fci, ^2, fcs) ■ 



(45) 



(46) 



If we use the leading-order solution ([271) for the mode function and neglect the variation of the terms C^'s for the 
integration of Eq. (|44p with the approximation a ~ — l/(iJr), the resulting bispectrum is [gH [6^ 



An D 



CiSi 



H 



4esF 2esF Mpi 8 4cf esF \Mpi 



2e,Fc2 V A/, 



H 
pi 



C7S7 



1 H 



(47) 



where 



Si 



S4 = 



5r 



i>j 



i>j 



(fclfc2fc3)' 
^3 



S3, St^^ ii + ^Y.'^^'^r 



3kik2k3 



1>J 



i>j 



2 

i>j 



Sh = ^ 



1 

if2 



- kik2k3 E - o fclfc2fc3if ' - 6 E - E + 9 E 



(48) 



^ Note that in Ref. [T^ the term 7?.^ is present in the expression of J^i , which gives rise to the first-order contribution r]s ■ We absorb this 
term to other coefficients, so that the field definition in Ref. [T3l is unnecessary. 
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and K = ki+k2 + k3. The five shape functions S'^'s (i = 1, • • • , 5) are present in the context of k-inflation [l^, [Sa]. In 
the Horndeski's theories the additional functions S'y and 5*8 appear, but they can be expressed by using other shape 
functions as [6^ 



5*7 — — -(35i — 5*2) 



18^3, 



— 3Si — S2 ~\~ 35*4 . 



(49) 



Since the three functions 5*6, 5*7, and 5*8 vanish in the hmit fcs — > 0, the last three terms in Eq. (j47p do not contribute 
to the local non-Gaussianities. 

The bispectrum A-jz coming from the contributions of Ci, C2, C3, Cg, C7 are 0-th order of e, while the bispectrum 
from C4 and Cs arc first order. Since the term C5 leads to the bispectrum at the order of e^, we can neglect its 
contribution. In the case where the leading-order terms oi A-ji vanish (which occurs for local non-Gaussianities), we 
need to take into account next-order corrections to the bispectrum coming from the integrals that involve the terms 
Ci, C2, C3, Cg, C7 in Eq. p2|) . Using the linear equation of motion SC2/STZ\i = 0, the Cg, C7, and Cs dependent terms 
can be absorbed into the first five terms in Eq. ([5^ [69| . Then the third-order action ([5^ reads 



53 = y dtd-^x^a-^CiM^iUn^ + aC2M^{R{&RY + a^C^MpilV + a:^ein{^^n){^^X) + a:\C^/M^y)d^n{dX) 
The coefficients Ci (i = 1, • • • ,5), which give rise to the corrections up to the order of e in A-r,, are [1^ 



C2 = C2 

C3 = C3 
C4 — C4 



(2-2e + ,.. + ,7-4.)C7-^C8 



H 



3c2Afpi 
3i? 



' -I- 3r/sF - 4s - 776)^6 + 



hi 



-(6 -I- 377sF - s - 2?77)C7 , 



C5 — C5 , 

where r]i = Ci/{HCi), with i = 6, 7, 8. The three-point correlation function analogous to (j47p is given by 



4e.^^ 



Ci5i 



1 



4e.^^ 



:C2>5'2 



1 ^ 



2e,F A/pi 



C3'5'3 + -C4S4 



(50) 

(51) 
(52) 
(53) 

(54) 
(55) 

(56) 



where we dropped the Cs-dependent term. The difference between Eqs. (|T7|) and ([55)) is that the bispectrum ([55)1 
includes the corrections coming from the time- variations of Cg and C7. However Eqs. (|47l) and (|56p are equivalent at 
leading order. The terms C,; (i = 1, • • • ,4) can be expressed as 

(57) 
(58) 
(59) 
(60) 







- ^ [e,<5Ci + (2e 


+ l^sF - 5?77)'5C7 + 3es(5C8] , 


C2 




- ^ [e,dC2 + (2e 
ci 


- risF - m + 4s)(5C7 - es^C8] , 


C3 


/rlcad 
— C3 i 






2 2 

3 3c^ 


C4 


2e,, 









where the leading-order terms are 





- --{ 


plead 




/rlead 
^3 


FMpi 

cm 



lead 



1 2A 
^ - 1 - TT 1 + 4(5C6 - -^SCj 



(61) 
(62) 
(63) 
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and SCi^s are the first-order slow- variation terms given by 

SCi ^ es-v,~4SG3x-l25Gix-S2SGixx + l2SG5^-10SG5x-SSG5xx, (64) 
SC2 = e,s + 77s - 2s + 4,5g4x + 2(5g5A' - 4(5g50 , (65) 

SC3 = —{Sg3X + 4(5g4X + 3Sg5X - ^40 - 4Jg50 + 8(5g4XX + 2Sg5Xx) - (3 + 2A3X) Sg3X -8(5 + 2\4x) Sg4XX 



- 4(4 + X5x)5g5XX + Sg44> + 8^50 - 8Sg4x - 9Sg5X - G— [(1 + A3x)(5g3X + ^(S^ 
5Ce = {I + )^3x)Sg3X + MS + '2X4x)SGixx + Sg5X + {5 + 2\5x)Sg5XX , 

= Sg3X + QSg4XX + dG5X + Sg5XX , 
SCs ~ Sg3X + 4:Sg4XX ■ 



(66) 

(67) 
(68) 
(69) 



In order to derive the full expression of Atz to the order of e, we need to compute the corrections to the first three 
integrations in Eq. ((50|). As studied in Ref. [HBl in the context of k- inflation, there are several corrections to the 
bispectrum ([55|. 

The first one comes from the variation of the coefficients Ci [i = I, 2, 3), i.e., 



C^{T) = C,{rK) - 111 — + • 

dt Uk tk 



(70) 



We evaluate all the physical variables at the time tk ~ —1/{Kcsk), which corresponds to the moment when the wave 
number K = ki + k2 + crosses the Hubble radius Kcsk = olrHk- 

The second one follows from the correction to the scale factor a ~ -~1/{Ht), i.e.. 



1 



+ \n{T/TK) + 0{e'). 



(71) 



Hkt Hkt Hkt 

Thirdly, the mode function ([TTI) is subject to change by taking into account the 0(e) terms on the r.h.s. of Eq. (|25p 

1 H y3/2 / I I 



u{y) = - 



l + -e+-s] e'^(^+5''=^)i7(i)[(l + e + 



(72) 



where y — Csk/{aH), v ~ 3/2 + e-|-7]si?/2 + s/2, and Hi^\x) is the Hankel function of the first kind. In the large-scale 
limit (y — > 0) the Hankel function behaves as Hj)'\x) — > — i/[sin(7ri')r(l — i>)]{x/2)^'^ and hence the mode function 
approaches 



.(0) 



1 



1 



2kV^ ^(esF),c,,Mp, fc3/2 



72 



72 



(73) 



where 72 = 71 — 2 -I- In 2 ~ —0.7296... and 71 = 0.5772... is the Euler-Mascheroni constant. Note that the quantities 
with the subscript k in Eq. (|73|) are evaluated at Cskk/ (ukHk) = 1. For the wave number ki there is the running from 
kj to K, as 



1 



ki 



1 



Hk 



1- e 



1 



VsF + 7;s] In 



Vi^sF)k, Csk, Mpi ^[esF)K CsK Mpi 
Writing the correction to the leading order solution ([27]) as Au*{t, k), it follows that 

1 



K 



Au*{t, ki) 



dT 



Au*{t, h) 



2k1'^^{esF)KCsK^^y 



(e -I- s){x — i) + ' 



+ 1 ( ^ + T^'^sF + ) (j — x) — isx" j In h y 2 



1 1 , 



3/2 



dH, 



(1)* 



diy 



(74) 



1^=3/2- 



,(75) 



2k^^^^{esF)KCsK^^P^ 
+1 e H — r]sF s + isx xm h 



xysx — le) + e + -?/sF + 77^ 



1 1 



d 

dx 



2 ' 2 ^ 

fdHl'-^ 



r3/2| 



\ dv 



=3/2^ 



, (76) 
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where x = —kiCsK t. 

In Appendix A we give the explicit forms of corrections to the first three terms in Eq. (|56p . Each correction can be 
expressed as 



4e,F 



lead 



K 



1 



lead 



4e,F 



K 



C'r"] SQs, (77) 



where 5Qi, 6Q2, and SQ3 are the 0(e) terms derived by summing up the contributions (|A1[) - (|A5|) . (|A6|) - (|A9|) . and 
(|XT0)) - (Pl3| . respectively. On using Eqs. (|6ip - (|63p . it follows that 



AAn = AA^^^ + AA)^' + AA-,^ 
1 f 1 



l(2) 



(3) 



7(3-1- ) {35Q, + 5Q,) + r(__i_ri:+4::^_ _in 1 



3/1 2A ^ JCe 4 SC^ 



4 V c' 



4 W 

The explicit forms of 3SQi + 6Q2 and (JQs are 
3SQ1+SQ2 

-2(7 + 271 + 672)6 + 3(1 - 271 - 2-f2)vsF - (21 - 1471 + 672)5 - 2(1 - 271)7)1 - 2(2e + rj^F + s) In 



(78) 
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i>j 



2(3 + 71 + 372)6 - 3(1 - 71 - ^2)VsF + ( 17 - y7i + 372 ]s + 2(1 - 71)771 



+(26 + TjsF + s) In 



kik2h 



± k^k^ + 3(271 - 1). E '^'4 + 1^ E - 3 i^'^)^ 



22e + 9r]sF + 2ryi + (11 + 471)5 



13 11 



kik2ks 



+ 7'?^^^ - 7'?i - 7^ E 



2 4 



i(3 + 71 +372)6- -^(1 -71 - j2)rhF + ^(1 -7i)?)i + 7(13 - 71 + 372)3 + ^{2e + T]sF + s)ln^^i^-^ 



^(1 + 7i)4 E - l^^'^ E ^^^1 + ^(26 + rj^p + .) (^2 E fcf + 2 5] fc,fc| - v) , 



(79) 



SQ3 



29 



2(2 + 271 + 372)6 + 3 - - 71 - 72 ri,F + 371 - 372 - — s + (271 - 3)7/3 - (26 + rj.p + s) In 



kik2kz 



(fclfc2fc3)^ 



e+\vsF + \A 



i^j i>j ' 



(80) 



where 77^ = {dCf'"^/dt)/{HCf'"^), and 



V = X - (/s| + fc^)Re 



U 



h{x) 



dxi ■ 



-iKxi / ki 



Xl 



- ^AA + perm. 

Kl 3 



-fciRc 



Xl 



dxi — ( fco + fc| + 7/02^; 



fc2 + fc a 



Xl e ''1 



t2+i3^,^dh*{xi) 



dxi 



{ki + fc3^)Re 



f 

Jo 



dxi- 



-iKxi I h\ 



X\ 



' 6 



r f°° 1 

> kf Rc / dxi 



fc2 + ''3 



''1 1 + 7- 



.^2+^3 



k^k 



Xl Tt"^! ) " (^'1) 



kl kl 



kl 

^Re 
kl 



^3/2 



dxi /i*(xi)e '-'1 



perm. , 



dxi h*{xi) 1 — * 



. k2 + ks 



kl 



Xl I e *=! " 



perm. 



dv 



'2ie^'^ + ie "(1 + 7x) [Ci(2x) + i Si(2x)] - iTrsinx + ittxcosx . 



(81) 
(82) 

(83) 



iy=3/2 



10 



The definition of Ai, G, and is given in Appendix A. In Eq. ((8T|) . we liave used tfie relation between the variables 
xj^ = —Kcskt and Xi = —kiCsKT = li2,3), as xk ~ {K/ki)xi (without summation over i). The variable 
xk = —Kcskt is related to Xi = —kiCsKT [i = 1,2,3), as xk = iK/ki)xi. The symbol "perm." stands for cyclic 
permutations with respect to fci, k2, and ^3. In Eq. (|79p we also used the relation 7)2 = ^1 + 2s to eliminate 772 (which 
follows from C^^*^ = — (Cs/3)C}°^'^). In Appendix B we evaluate the values of V and U as functions of r2 = fc2/fci and 
rs EE ks/ki. 

The total bispectrum A-r, is the sum of Eqs. (|5S)) and ([75]). which can be written as 



lead , yi corre 



where 



A lead 



A corre 



4 
1 
4c? 



1 SCj 



(35i - 52) + 
(5Ci + (2e + 7?7,j. - ^m)— + 3(5C8 



3/1 A 3A ^ 6(5C6 



6 



(84) 
(85) 



(5C7 



(5C2 + (2e - ?7,F - 7^7 + 4s)— - <5C8 



So 



SC. 



1 



SC7 



-6C3 + {3t]sF - 7?6 - 4s) ni^VsF - 2777 - s) 



5*3 - -r^i'^s - 3(5C8)S'4 
4c„ 



1 ( 1 



2 JC7 



1 



2A 



4 (5C7 



(86) 



The leading-order bispectrum A^^^ (given already in Refs. [61I, [g^I) and the correction AS^^'^ are of the orders of 
0(6°) and ©(e), respectively. 

IV. LOCAL, EQUILATERAL, AND ENFOLDED NON-GAUSSIANITIES 

The non-linear parameter characterizing the strength of non-Gaussianities is defined by 

10 An 



NL 



3 Ell n 



(87) 



In the following we estimate /nl for three different shapes of non-Gaussianities. 



Local non-Gaussianities 



The local shape corresponds to ^3 — > and ^2 — > fci = fc, in which case /jl^L^' = [^jZ^A-nlk? . Since 5*1 = ^"^72, 
S-2, = 3A;^/2 = 351, and 5*3 = 5*4 = 0, the leading-order bispectrum ([551) vanishes. In the limit that fc3 — > the 
function U given by Eq. approaches fc'^/2 [H^ (see also Appendix B), so that the term (5(53 in Eq. ((5(1)) vanishes. 
Then the bispectrum (|86p reduces 



A^n"" = |^('5Ci + 3,5C2) - - 1) (35^1 + 5^2) + ^ [4(2e + 7;,^^ - 2r;7 + 3s)fc3 + 4(3,5Qi + 



Using Eqs. ([M)) and (|55|) together with the relations p^ . (pTT). and (^5)) . we have (5Ci -I- 3(5C2 = 4e + 277^^? — 6s. In 
the limit fc3 — > the function V behaves as V — > 20A:'^/3 [5^1 (see also Appendix B), so that Eq. ([79| reduces to 
3(5gi + (5Q2 ^ (2e - ?7sF + 5s + 2f]{)k^ . Then Eq. (HH]) reads 



A^ 



corre 

n 



fc3 



_fc3 
25" 



(2e - 77^F + 5s 2^71) + T^(77sF - 4s - 771) + k^^{2e + As + fji- 777) . 



(89) 



Taking the time-derivatives of Ci and C7, we obtain the following relation 



^Cj 1 1 

^— (2e + 4s + 771 - r/7) = -{i]sF 2s) - -T^ijlsF - 4s - ^7i) 



(90) 



Substituting Eq. ([90]) into Eq. ([89]), it follows that 



^™- = _(2e + 77.f +s). 



(91) 
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Using the spectral index n-ji given in Eq. (|29p . the non-linear parameter is expressed as 



./nl"' 



"°^^' = T^t(1-"^)- (92) 



12 

This matches with the Maldacena's result [l^ derived for a canonical scalar field (see also Refs. [s^ |H, for the 
derivation of the same relation in other single field models). Creminelli and Zaldarriaga [13] pointed out that the 
consistency relation (|92l) should hold for any slow- variation single-field inflation. In fact we have shown that this holds 
for most general single-field scalar-tensor theories with second-order equations of motion by explicitly computing the 
slow- variation corrections to the bispectrum ([55| . Since I/^l'^'I is much smaller than 1 in such models, the observational 
detection of local non-Gaussianities with I/^l^'I ^ 1 implies that we need to go beyond the slow- variation single-field 
scenario. 

In the limit that ^ Q the shape functions Sg, Sj, and vanish. This means that the functions Cg, Cj, Cg on 
the r.h.s. of Eqs. ([?T|) - ([5ll) do not contribute to the local non-Gaussianities. In fact, we can derive the consistency 
relation ((92|) by setting 5C^ = in Eqs. (|88p and (|90)) . In this sense the situation is analogous to that in k- inflation. 

Let us consider the not-so squeezed case in which the ratio = k^/ki is no n- vanishing, i.e., < ra ^ 1 and 
ki ~ k2. The leading-order non-linear parameter following from Eq. (j85|) is given by 

.lead _ 5rg [ai(22 + 4r3 - 4r§ - ri)(2 + r,) + 2a2] 

/nL — o^o , „ \Rfn , „3^ ' y^'^) 



where 



3(2 + 7-3)3(2 + ri) 
1 1\ ^ 1 SCr 3 f 1 ^\ 3A , 65Ce 6 SCj 



In the regime ra ^ 1 we have /nl'^ — 5(22q;i + Q;2)r|/24. The values of ai and a2 depend on the models, but for 
Cg ^ 1 they are at most of the order of 1/c^. In this case the leading-order non-linear parameter can be estimated as 

Then the transition from the value (1921) to the value (l95l) occurs at 



ra w Cs VI - n-jz . (96) 



The effect of the term l/c^ in Eq. becomes important for > Cg y/1 — n-ji. When n-ji ~ 0.96 this condition 
translates into > 0.2 Cg. If Cg = 0.1 and > 0.1, for example, the non-linear parameter (|95p can be larger than 



the order of 1. However, for Cg = 0.1, we also expect the dominant contribution to /nl to come from other shapes 
(equilateral, orthogonal, etc.) 

For the models in which is close to 1, ai and 012 are of the order of unity. Hence the leading-order non- linear 
parameter can be estimated as ~ t"! <C 1 in the regime ra ^ 1. By increasing the value of ra from 0, we can 

observationally discriminate between the models with 1 and w 1. 



B. Equilateral non-Gaussianities 

The equilateral shape is characterized hy ki = k2 = k^ = k, in which case the non- linear parameter is /^l"'' = 
{10/9){ATz/k^)- Since = 4k^/3, S2 = 17k^/6, and S3 = k^/27, the bispectrum ^ gives the leading-order 
non-linear parameter 

_^equil,lcad ^ g _ _ i^^ + [(1 + X,x)Sg3X + 4(3 + 2Xix)SGiXX + Sg5X + (5 + 2X5x)Sg5Xx] 

65 

{Sg3X + (3Sg4XX + Sg5X + Sg5Xx) ■ (97) 



162c2e 



If c2 « 1, then we have » 1 



12 



In the equilateral limit the functions V and U are given by V = 15[1 + (1/2) ln(2/3)]fc3 and U = [61n(3/2) - l]fc3, 
respectively (see Appendix B). Then the functions 5Qi2 = (3i5(5i + 5Q2)/k^ and 5Q^ = SQ^/k^ reduce to 

SQ12 = e[6-781n(2/3)-147i]/3-r;,F[77i + 261n(2/3)]/2 

+s[16-391n(2/3)]/3 + ?7i(77i + 20)/6, (98) 
SQ3 -e[107i + 55 + 1681n(2/3)]/27-7/sF[27i +7 + 281n(2/3)]/9 

-4s[10 + 21 ln(2/3)]/27 + 7)3(271 - 3)/27. (99) 

The correction to /^'l"' '"^'^ coming from Eq. ((89|) is given by 



equil.corrc ^ 

~ 972^ 



3e^(24(5Ci + 5ISC2 + 4c,(5C3) + 8c^(3r;,,_F - 4s - i]e)SCe + (450e + 327?73i=^ + 620s - mr]j)6C7 



For the theories in which /^i"'''''^'''^ vanishes, the next-order correction /^'l"'''^""'' is the dominant contribution. In the 
case of a canonical scalar field with the Lagrangian P = X — V {(/)), G3 — 0, G4 — 0, G5 = 0, for example, it follows 
that = and = 55e^^^^ 

For the theories with 7^ 1 the non-linear parameters ([97|) and (jlOOp reproduce the results known in literature for 
sgecific models of inflation. For example, this is the case for k- inflation [ijj HI] : k- inflation with the Galileon terms 



62 , |63[ , potential-driven Galileon inflation |73l | , and inflation with a field derivative coupling to the Einstein tensor 
53|. Generally we require that ^ 1 to realize the large equilateral non- linear parameter. 



C. Enfolded non-Gaussianities 



The enfolded shape is characterized by ^2 + ^3 — ki. Taking the momenta fci = fc and k2 — )■ k^ — fc/2, the 
non-linear parameter^ is = SAn/i^k^)- Since = 23k^/64:, S2 = 63fcV64, and 5*3 = k^/128 in this case, the 

leading-order non-linear parameter is given by 



^cnfoldjcad ^ ^ ^ M _ + _L [(1 + X,x)Sg3X + 4(3 + 2X,x)SGiXX + W + (5 + 2X,x)Sg5Xx] , (101) 



where, unlike the equilateral case, the (SC7-depcndent term in Eq. (j85p disappears. 

In the enfolded limit one has V = [315/64- ln(2)/2]fc3 and U = [5/24 + ln(2)/4]fc^ (see Appendix B). The functions 
6Qi2 and SQ3 are 

SQ12 = e[45 - 2471 + 108 ln(2)]/64- ry,,F[87 + 3671 - 108 ln(2)]/128 

+s[333 + 1081n(2)]/128 + r)i(39 + 37i)/32 , (102) 
SQ3 = -e[157i + 17-541n(2)]/192-?7,_F[367i-5-1081n(2)]/768 

-s[109 - 1081n(2)]/768 + %(27i - 3)/128 . (103) 

ine correction to is 

Ilfold f^rw-re' 1 



^entom.corre ^ ^_[g^(23(5Ci + 63SC2 + SciSCs) - 2ci{m + 4s - 3rjsF)SCe + (172e + 254s + 92r;,i. - 17irj7)6C7 



For a canonical scalar field we have that j^^f°''^''<=^'^ = g and j^'^f°i'i>'^°"'<= — 7es/8 + 577s/12. 



At the point k\ = k and fc2 — > fea = fc/2, the equilateral shape gives no contribution. However, the contribution from the orthogonal 
shape will be of the same order of the enfolded one, since, by definition, in this case we have ^ ^2/^'J^'°''^. 
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V. SHAPES OF NON-GAUSSIANITIES 



The leading-order bispcctrum (j47p can be written in terms of the sum of each component, as A-jz — X]i=i ^tz 
Then the function J--ji{ki, k2, fcs) in Eq. (j46|) is decomposed into eight components 



n 



A. 



n 



(105) 



where BiS are the coefhcients appearing in front of each shape function Si in Eq. (|T7)) . say, Bi = c^CiJjAesF). 
In order to estimate the correlation between two different shapes, we define the following quantity |42| 



0) TrW)^ 



(106) 



where 



k2, h) 



(107) 



The integration should be done in the region < fci < oo. < k2/ki < 1, and 1 — ^2/^1 < k^/ki < 1. Note that 
the above integral can be expressed in terms of r2 = k2/ki and r^, = k^/ki with the integral of fci factorized out. 



For \C{J^!^ ,J^!ji^)\ close to 1 the correlation is large, whereas for \C {J-!jl' , J-!jl' )\ close to the two shapes are almost 
orthogonal with a small correlation. 

The CMB data analysis of non-Gaussianities has been carried out by using the factorizable shape functions which 
are written as the sums of monomials of ki, k2, and k^ . There are a number of templates J^-jz which resemble model 
predictions of the bispectrum. The templates corresponding to local and equilateral non-Gaussianities are given, 
respectively, by [13, EH 



J-]^=^'(fci,fe,A;3) = (27r) 



(108) 



and 



^r\kuk2,k,) = i2n)^{^f^r' 



;,3l3 i.31,3 i,3 2,3 l,2jL2L2 I L l2 l3 
1 2 2 3 3 1 <v-|^A(2'i'2 \ l 2 3 



5 perm. 



(109) 



Since the local non-Gaussianities are small in the Horndeski's theories, we do not consider the correlation with the 
local template. 

The orthogonal template, which has a small correlation with the equilateral one, is given by (isj 



770rtho 



(A:i,fc2,fc3) = (27r)4 



1,31,3 i,3l3 l3i,3 1,21,21,2 ' I 1,, 1,21,3 

12 2 3 3 1 iv-j^ r\i'2 iv^ \^ rv]^ rv^ hj^ 



5 perm. 



(110) 



The enfolded template, which is a linear combination of the orthogonal and equilateral templates, is defined by [44 



J-^*"^(fcl,fc2,fc3) = (2^)^(^/^lf 



1 



1 



1 



p)L3 
"'1"'2 



1,3^,3 
"'2 "-3 



L3p ' Ju2t,2jL2 
3 1 12 3 



1 



^1 ^2 ^3 



5 perm. 



(Ill) 



In Table U we show the correlation between (i = 3,4,5,7,8) and the three templates given above. Since 



the correlations between t!)^ (i = 4, 5, 7, 



and are close to 1, {i 

the equilateral shape. In particular, both C and C(J"4*\ 7"^''*'^°) 
correlation between J^^"'' and 



- 4, 5, 7, 8) are well approximated by 
are close to with the same level of 
Hence the shape functions 6*7 and 5*8 highly mimic the equilateral template. 
They also vanish in the local limit (fcs — > 0) and in the enfolded limit (^2 + ^3 — > k\). Note that both \C(T^ 
and |C(J"^^\ J"^'^f°''i)| are close to 0.5 (which is similar to the value \C{J^^'^'\T^^°^% = 0.511911). Since there is the 
relation = (T"^"^'^ - F^^^°)/2, we have |C(J"W^ j-«ifoid-)| ^ o.5 for the shape function f'^ very close to the 



equilateral one. 
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T^cquil 
■^TZ 


X*ortlio 


X"cnfold 


■^ TZ 


TZ 


■^ TZ 


T-(7cquil) 
■^TZ 


■* TZ 


T-(7ortho) 
■'TZ 


oTr(i) 7r(2) 
^■'TZ ■'TZ 


•^^quii 


1 


0.0254062 


0.511911 


0.936177 


-0.998757 


-0.994234 


0.999892 


-0.999994 


-0.00693357 


0.986954 


-r*orttio 




1 


-0.845755 


-0.290742 


0.0177139 


-0.117961 


0.0353534 


-0.0277283 


0.904843 


-0.116557 


-r-cntold 
TZ 






1 


0.749518 


-0.548302 


-0.4293 


0.503306 


-0.509913 


-0.781246 


0.626939 


71 








1 


-0.952469 


-0.893224 


0.933797 


-0.935615 


-0.357802 


0.980504 












1 


0.987653 


-0.998384 


0.998686 


0.056402 


-0.993745 














1 


-0.994696 


0.99437 


-0.0995524 


-0.964012 


T-(7cquil) 

■'u 














1 


-0.999936 





0.9859 


■'TZ 
















1 


0.00524865 


-0.986715 


-^(7ortho ) 
■''TZ 


















1 


-0.167335 



Table I: The correlation (|106|l between two different shape functions. J^.j^'^''"''' is the normalized shape of , whereas J-] 
is the shape function orthogonal to JT^'^''"''' (or, equivalently, to -^-J^'). 



We recall that the functions 3Si — S2 and 5*3 are related with 5*7 [see Eq. (|49|) ]. We can use this property in order 
to rewrite the leading-order part of the bispectrum in a convenient basis. We introduce the following shape 



1 9 

rrcquil c 



(112) 



whose minus sign has been chosen so that has a positive high correlation with the equilateral profile. Fur- 

thermore, we can analytically show that the following shape is exactly orthogonal to 

12 



■cquil _ 



ortho 



14 - 13/3 



(^57 + 3Si - ^2) 



where 



/3 



16 248041 - 252007r2 
y 1986713 - 2016007r2 



1.1967996. 



(113) 



(114) 



The normalizations of ^y^"'' and S'y'"^" have been done such that, at the equilateral configuration (fci = fc2 = ^"3 = k), 
we have 5^''"'' = = 5*7'*''°. This normalization follows from the standard definition of the previous templates 
introduced in the literature. 

We note here that the leading-order bispectrum (|85|) includes the term 35*1 — S'2, so that we also consider the 



correlation between the combination 



(1) 

TZ 



correlation with Compared to J^!^'' and J^!j^' , however, it is not very close to the equilateral shape. Moreover 

has some correlation with the orthogonal shape, i.e., C(J^^'', J^^'*''°) = —0.290742. 



J-".^'' and other shapes in Table L The shape SJ-"^"* 

-(8) 



t!^^ has a high 



As we see in Table HI the correlation between J^j^ 
functions 35i — S'2 and ^3 are expressed as 



(7ortho) J".^"'' is very small. In terms of S'^''"'' and £'7''"^°, the 



35*1 — 5*2 



Q rrcquil 

12^*^ 



14 
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(115) 



Using these relations, the leading-order bispectrum 
the orthogonal basis 3^'-'^^°, as 



432 ' ' ' ' 144 
can be written in terms of the equilateral basis S'^'^"'' and 



Alcad _ qcq 



C25? 
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where 
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C2 



13 
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24 



(2- 



■3^) + Y^(2-3/3) 



SCe 

6ec 



(2 - 3/3) 



3e^c2 



14- 13/3 
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A 

is 
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(116) 

(117) 

(118) 



The coefficients ci and C2 characterize the magnitudes of the three-point correlation function coming from equilateral 
and orthogonal contributions, respectively. Finally, we also introduce the enfolded shape function 



enfold 



(^7^1"' -5?'-*^°)/2 



(119) 
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Figure 1: The shape functions J"4^"''""'(l, fca/fci, fc3Ai)(fc2/fci)^(fc3/fci)^ (top left), T^^°"''^°\l,k2/ki,ki/k^){k2/k^f{k3./k^f 
(top right), and T!j^'^''^^°^'^\l, k2/ki, k3/ki){k2/ki)^ (ks/ki)^ (bottom). 



which has a maximum at ki 2fc2, ^"2 ^3- Note that 5'™'^°''^ vanishes at the equilateral configuration. 

In Fig. [1] we plot the three shape functions and jrO'*'"^"^'^) multiplied by the functions 

{k2/kiY{h/kiY. The correlations of P^-"^ = (27r)^yl^^'^/ JlLi with the shape functions -F^'*"'', and 

■cnf 

n 



jronfoid ^^-^ evaluated as 



., , , 8.25104 (1.44717 X lO-^ci- 1.70045 X 10-%) 

^ ^1.42610 X 10-2c2 + 4.09480 x IQ-^c^ 

, , . 5.12494 (8.23792 X 10-% +3.57273 X 10-2c2) , ^ 

^1.42610 X 10-2c2 + 4.09480 x IQ-^c^ 
^enfold ^ ^._^lead _^c„fold^ ^ 8.80788 (6.82395 X lO'^d - 1.79487 x 10-2c2) 

^1.42610 X 10-2c2 +4.09480 x lO^^cS 

which depend on the coefficients ci and C2. In particular, we find C(J-'j^'^"^°''^'', = 0.928621. 



VI. SHAPES OF NON-GAUSSIANITIES IN CONCRETE MODELS 



Let us study the non-Gaussianities of concrete models of inflation in which the bispectrum (|116p can be large due 
to the small scalar propagation speed Cg. As we will see below, there are some models where the orthogonal shape 
provides an important contribution to the bispectrum. 
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A. Power-law k-inflation 



We first consider k-inflation characterized by 

Pi(b,X) = K{<j)){-X + X^), G3=0, G4-O, Gs^O, 
where K{(t)) is a function in terms of (p. From the background equations ([9|) and ([TOl) it follows that 



(123) 



3Af2jii-2 
X{3X - 1) 



i/ _ 3(2X - 1) 
- 3X - 1 



(124) 



As an example, we study power-law inflation characterized by a oc t^^'^ and H = l/{'-ft), where 7 (^ 1) is constant. 
Substituting the Hubble parameter into the second of Eq. (|124p . we obtain X ~ {3 — 7)/[3(2 — 7)] and (f) = (po + 
^(3 - 7)/[3(2 - 7)]t (00 is the initial value of the field). From the first of Eq. ([T^ we find that power-law inflation 
is realized for the choice 1541 



6(2 - 7)Mgi 



In this case and A/S in Eq. (|116p are given by 



7 



3(4 - 7) 



A 1, , 



For 7 < 1 we have that l/c^ ~ I2/7 > 1. 

In this model (dubbed Model A) the leading-order bispectrum (|116p reduces to 



(125) 



(126) 
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(127) 



In the hmit < 1 it follows that A^^'''^ ~ {-0.252/cl)S^ 



-fcquil 



(0.016/c2)5? 



■tho 



parameters ((97)) and pOip are given by /^i"'''''^'''^ — ~85/(324c3) and ■f^°^'^'^'^^'^ r- 
WMAP9 year constraint ^ of the equilateral shape gives the bound > 1.2 x lO"-^ (95%CL). 
The scalar spectral index (|29p and the tensor-to-scalar ratio (pij) can be expressed in terms of Cg 



Since in this limit the non-linear 
l/(32cs) respectively, the 



nn = l 



24cg 
l + 3cr 



192c3 
l + 3c? 



(128) 



For the ACDM model without the running scalar spectral index, the bounds on Ug and r from the WMAP9 data 
alone are nn = 0.992 ± 0.019 (68 % CL) and r < 0.38 (95 % CL), respectively. If we combine the WMAP9 data with 
the measurements of high-^ CMB anisotropics, baryon acoustic oscillations, and the Hubble constant, the constraints 
are 0.9636 ± 0.0084 (68%CL) and r < 0.13 (95%CL). If we employ the bound 0.95 < < 1 then the scalar 
propagation is constrained to be < 2.1 x 10""^. Since in this case r < 0.018 from Eq. (|128p . the observational 
constraint on r is satisfied. Combining the bound of cj. with that of the scalar non-Gaussianity, it follows that 
1.2 X 10-3 < c2 < 2.1 X 10-3. 

In Fig.HJwe plot the shape function ^"^^'''^(1, k2/ki,k3/ki){k2/ki)^{k3/ki)^ for = 2.0x10"^ (labelled as "A" in the 
figure), where J"]^^^ = (27:)^ A^^""^ / OLi ^f- When = 2.0 x lO'^, the correlations (fT^ - p^ are G^^"" = -0.99474, 



ortho 



0.06305, and C 



enfold 



-0.58511, respectively. Hence the shape of non-Gaussianities is close to the equilateral 



one illustrated in Fig. [TJ whose property is is independent of the choice of 1 



B. k-inflation with the term G3{X) 
We study k-inflation in the presence of the covariant Galileon term G3{X) characterized by [63l. iTlj 



P{X) = -X 



2Af4 ' •^31,-^; ^4, 



(129) 
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where M and /x are constants having a dimension of mass. In this model (dubbed Model B) the de Sitter solution is 
present when e = 5px + ^^g^x = 0- Using Eq. as well, we obtain 



18/l2 X 



(130) 



where x = X/M^. As long as inflation is realized in the regime n/M-p\ ^ 1, x is close to 1. In what follows we replace 
X for 1 except for the terms including 1 — x. 

Along the de Sitter solution we have that A/S = Gc^/e,,, 6Cq = 5C7 = (5g3X = 2(1 — x). Since = 2(1 — x) and 



(1 — x)/6, the bispectrum (|116p reads 
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ortho 



(131) 



If < 1, we obtain Ai^ 

2\ ] ^enfoldjead 



lead 



/•cquil,lcad 



(0.109/c2)S'7'i"" + (O.Oie/cf)^^ together with the non-linear parameters /^l 
5/(36c2) and /«'^'°^'^''«'^^ ~ -l/(32c2). Notice that the sign of /nl"' ' is opposite compared to the k-inflation model 
A. From the WMAP 9 year bound ^ it follows that > 4.3 x 10"^ (95 % CL). In Model B the scalar spectral index 
is n-ji = 1, whereas r is related to via 



3V3 
48 ' 



.2/3 



(132) 



The WMAP9 bound r < 0.38 gives the constraint < 1.6 x 10 ^. Using the severer bound r < 0.13, it follows that 
cl < 7.7 X 10""^. In both cases there are viable parameter spaces compatible with the constraint from the scalar non- 
Gaussainity. When = 2.0 x IQ-^, for example, we have C^^"" = 0.96865, C"**^" = 0.25272, and C™f°''^ = 0.29984. 
As we see in Fig. [21 the shape of non-Gaussianities for = 2.0 x 10"^ is approximately close to the equilateral one. 

For potential-driven inflation {P — X — V{(f))) with the term G3{X) = iiX/M'^ [t^I, the non-Gaussianities are small 
because is not much smaller than 1 fTSj . In such models, if the Galileon self-interaction dominates over the standard 
kinetic term even after inflation, there is an instability associated with the appearance of the negative after the 



field velocity cj) changes its sign [7J] . In the model B discussed above, reheating needs to occur gravitationally [71 
so the situation should be different from that in Galileon inflation driven by a potential with a minimum. 



C. k-inflation with the term G4{X) 



The next model (dubbed Model C) is k-inflation with the covariant Galileon term 6*4 (X) [64I, i.e. 



A2 

2M^ ' ^^^^> " 
Similar to the case (jl29p . there is a de Sitter solution satisfying the conditions 

1-x fiM _ 1~ X 

'mJi ~ 6a;2(3 - 2x) 



P{X) = -X 



(133) 



(134) 



36// X 

Inflation occurs in the regime where x = X/M"^ is close to 1. Employing the similar approximation to that used 



previously, we have SCe = 2SCy = 12(5g4xx, £s 
Hence the bispectrum (|116|) reduces to 



AXX, OgaXX 
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ortho 



(135) 



{0. 018 /c^,)Sj'^'"^ + (0.016/c2)S'f The equilateral and orthogonal non-linear 
parameters in this limit are given by /^l"'' ''^^'^ — 25/(648c^) and ~ — l/(32cs), respectively. The tighter 



In the limit ^ 1, we have A^-^'^ 



^ In the framework of the effective field theory the Lagrangian is valid at the energy scale of inflation, but some higher dimensional 
operators can appear during the reheating stage. There may be a possibility that the Laplacian instability can be avoided by such 
operators. 
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Figure 2: The leading-order shape function T^'"^{l,k2/ki,k3/ki){k2/ki)^{k3/ki)^ for = 2.0 x 10 ^ normahzed to 1 at the 
equilateral configuration. Each panel corresponds to the models A, B, C, D. 



constraint on comes from the enfolded bound rather than the equilateral bound i.e., > 3.6 x 10~^. Note 
that the scalar spectral index is n-ji ~ 1 and that the same relation as Eq. (|132p holds in Model C. Hence the upper 
bound of ct coming from the observational constraint of r is the same as that of Model B, i.e., < 1.6 x 10"^ for 
r < 0.38 and < 7.7 x lO'^ for r = 0.13. 

In Fig.[2]we plot the shape of non-Gaussianities for = 2.0 x 10^-^, in which case the shape is between the equilateral 
and orthogonal ones shown in Fig.Hl In fact, we have C"=i"'i = 0.58145, C°''*^° = 0.75322, and C™f°''i = -0.33691 for 
Cg = 2.0 x lO^'^. The orthogonal contribution tends to be less important for the values of larger than the order of 
10-2. 



D. k-inflation with the term G5(X) 
Finally we study the following model (dubbed Model D) 



P{X) = -X 



2M4 



G5{X) 



(136) 



In this case, for x — X/M^ close to 1, there is a de Sitter solution satisfying the conditions = Af^/(6il/p[) and 



/i^MVAfS = 27(l-a;)V25. Since 5Ce = 3(5C7 = 36(1 -x)/5, = 18(1 -x)/5, cl = 3(l-a;)/10, and A/S = 1/2, the 
bispectrum is given by 
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Figure 3: The correlations (fT20ll - (fT22t versus for Model D. In the limit that < 1, IC""'"'"*! and |C°'th°| are larger than 
|C""J""|. On the other hand, for cl > 0.04, |C"='»""l is largest. 

In the limit c2 < 1 we have A^^'^ ~ (-0.012/c2)S'7''"" + {OmQ/cl)Sf^'^° = -(0.023/c2)S'f"f°i'^ + (O.OOS/c^) 5'f *ho^ 
in which case the sign in front of the shape S'y'^"'' in Eq. (|137p is opposite to those in the models B and C. In the 
same limit the equilateral and enfolded non-linear parameters are /^l"'' ''^''"^ — 5/(9720^) and ~ — l/(32c^), 

respectively. In this case y^'i^"'^''®^'* jg smaller than t>y one order of magnitude. The WMAP 9 year enfolded 

bound (U) gives the constraint > 3.6 x 10~^. In Model D we have that n-ji = 1 and that the relation between r and 

is the same as Eq. (|132p . Hence the upper bound on is the same as that of Models B and C. 

The shape of non-Gaussianities for ~ 2.0 x lO^'' is plotted in Fig. [2] When = 2.0 x 10"'^ the correlations 
(fT20| -(fT22 l) arc C"=^"" = -0.35587, C°''"^° = 0.83545, and 6"="^°'^ = -0.90787, respectively, in which case the shape 
has quite high (anti)-correlations with both the orthogonal and enfolded templates. In Fig. [3] we show the correlations 
(-.cquii^ C°'*''°, and (7™^°''^ versus c^. For > 0.04 the correlation with the equilateral template is larger than those 
with other templates. For < 0.04 the contributions of the orthogonal and enfolded shapes tend to be important. 
In the limit that <C 1, |C'°if°i'^| is largest among other correlations. Model D is an explicit example where the 
orthogonal (or enfolded) shape provides a significant contribution to the bispectrum. 

VII. CONCLUSIONS 

In the Horndeski's most general scalar-tensor theories we derived the three-point correlation function of primordial 
curvature perturbations generated during inflation in the presence of slow- variation corrections to the leading-order 
bispectrum. Unlike previous works [6ll l62j. the bispectrum (j84p is valid for any shape of non-Gaussianities at first 
order of e. 

In the squeezed limit [k^ — >■ 0, fci — )■ ^2) the leading-order bispectrum (1851) vanishes, so that the correction ((86)) is 
the dominant contribution to A-r,. By using Eq. (|86p. we showed that the non-linear parameter in this limit is given 
by /nl'^' = (5/12) (1 — n-ji). This agrees with the result of Refs. [HI, |43| in which the three-point correlation function 
was derived by dealing with the long-wavelength curvature perturbation (mode ^3) as a classical background. As 
demonstrated in Ref. [47| , this result should be valid for any single-field infiation in which the decaying mode of TZ is 
neglected relative to the growing mode. Our direct computation of the three-point correlation function in the presence 
of all possible slow-variation corrections is another independent proof that the non-Gaussianity consistency relation 
holds for most general single-field theories with second-order equations of motion. 

The result of the local non-Gaussianities shows that I/^l^^I is much smaller than 1, e.g., f^ff' = 0.0125 for 
n-ji = 0.97. In the case where the shape of non-Gaussianities is not exactly the squeezed one (0 < rs = fcs/fci <g; 1), 
the leading-order bispectrum gives the non- linear parameter \f^ff\ ~ rl/c^. Hence the leading-order term dominates 
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over the correction for > Cgy/l — n-ji. For the models with ^ 1 the non-Hnear parameter can be as large as 
|/nl| > 1 with the growth of r^. By measuring the shape which is not so squeezed, it should be possible to discriminate 
the models with different values of cj.. The leading-order non- linear parameters in the equilateral and enfolded limits 
are given by Eqs. ([W)) and (|10ip . respectively, whose magnitudes can be larger than 1 for the models with <C 1. 
These results will be useful to constrain concrete models of inflation in future high-precision observations. 

We also showed that the leading-order bispectrum can be expressed in terms of the sum of the two bases S'^'^"'' and 
^ortho^ The shape S'^'^"'' is very highly correlated with the equilateral template (|109|) . It also vanishes in both local 
and enfolded limits. The shape 3™*"^°, which is defined by (|113p . is exactly orthogonal to S'y'^"''. The coefficients ci 
and C2 in front of S'^'^"'' and S™^^^° in Eq. (|116p characterize the equilateral and orthogonal contributions, respectively. 

In Sec. I VII we presented concrete models in which the orthogonal shape can provide important contributions to 
the bispectrum. In power-law k-inflation the shape of non-Gaussianities is well approximated by the equilateral 
type. However, in k-inflation described by the Lagrangian P{X) ~ —X + X'^/{2M'^) with a number of different 
Galileon terms like G4{X) = /xX^/M^ and G5{X) = /xX^/M^°, we found that the orthogonal contribution is crucially 
important for ^ 1. In the presence of the term G5{X) = iiX'^/M-^'^, the correlations with the orthogonal and 
enfolded templates in the regime <C 1 are higher than that with the equilateral template. 

It will be interesting to see how the observations such as Planck [t^ provide the constraints on the scalar non- 
Gaussianities as well as the scalar spectral index and the tensor-to-scalar ratio. In particular, if future observations 
confirm the value > 1 at more than 95 % CL, this implies that we need to go beyond the slow- variation single- 

field inflationary scenario (including the Horndeski's theories). The information of other shapes of non-Gaussianities 
(including the not so squeezed one) will be useful to discriminate between many different models. 
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Appendix A: Corrections to the bispectrum 



Following the calculations of Ref . [5al we present the explicit forms of the corrections to the bispectrum ([55]) . They 
come from the first three integrals in the action ([50)) . We write each OjV) contribution to the bispectrum as A^^"* 
(i = 1, 2, 3). Note that similar calculations were also carried out in Ref. j3l| . 

(i) / dt(fxa^CiM^{Rn^ 
• (a) The correction from the variation of Ci in Eq. (|70p 
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(b) The correction from the scale factor a in Eq. ([TT]) 
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(d) The contribution from the correction to u*{T,ki) in Eq. ((75l 
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(e) The contribution from the correction to -^u*{t, ki) in Eq. (|76p 
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where = ^Kc^j^t, and 
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(b) The correction from the scale factor a 
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(c) The contribution from the correction to u(0, fcj) 
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(d) The contribution from the correction to u* (r, ki 
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where is defined in Eq. 



Appendix B: Expressions for U and V 

We evaluate the integrals U and V defined in Eqs. ([5^ and ((ST|) for general values of r2 = fc2/^"i, '"3 = k^/ki, and in 
three different limits: local shapes {r^ — 0, r2 = 1), equilateral (r3 = l,r2 = 1), and enfolded {r^ 1/2, r-z — > 1/2) 
shapes. In doing so, we need to use the following relations 

h*{x) = 2sin(x) + [sin(a;) — xcos(a::)]Ci(2a:) — [cos(x) + xsin(a::)]Si(2a;) 

+i {2 cos(x) + TT sin(x) — nx cos(x) — [cos(x) + x sin(x)]Ci(2x) — [sin(a;) — x cos(x)]Si(2a;)} , (Bl) 

dh*(x) 

— = cos(x) — sin(x)/a; + a;sin(a;)Ci(2a;) — a;cos(x)Si(2a;) 

dx 

+i [7ra;sin(x) — sin(x) — cos(a;)/a; — a; cos(x)Ci(2a;) — a; sin(a;)Si(2x)] . (B2) 
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The U integral 



Let us first evaluate the integral U. We employ a few different procedures of regularizations, but they lead to the 
same final results. One possibility is to solve directly the integral, as 
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where fu is the integrand in Eq. ((82|) . We can either neglect the terms which rapidly oscillate around 0, or we can, 
equivalently, shift y in the complex domain for the oscillating terms like e^**^^ — >■ e^'*^'' (i±'/p) q (Method I). 
Another possible method (Method II) consists of solving the following limit 
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We need to check whether both methods lead to the same results, and we can - at least for some values of r2, this 
proves to be possible - check the result against numerical integrations. 

By employing the Method I mentioned above, the result of the integration is 



2rlrl 



In 2 



Arlrl (3r| - 6r2 + r| + 3) 
(r2-r3-l)3(r2+r3-l)3 " (ra-rg + l)^ 

^1 / [(^2 + 1)^ - 6r§ (r2 + 1) + rf]r| r|[rf + Srgr^ + 3(r| - 2)r2 + r| - Grg + l]r| 



{r2 + ?-3 + 1) H (»-2 - 7-3 + 1)2 (r2 + r3 - 1)2 

rl[rl - 6 (r3 + 1) rg + (r3 + 1)^] ^ 2r2 (ra + ra + l)2r2 ^frj^r^+V 

(1 - 7-2 + r^) 2 (r2 - r3 + 1) 3 \^ r3 

2r2 (r2 + 7-3 + 1) 2r2 ^r2 + r3 + 1 ^^ ^ 2r2 (7^2 + r3 + 1) 2r2 



(1 - r2 +7'3)- 



r2 



{r2 + r3 - 1) 3 



In(r2+r3 + l) 



(B5) 



where U{ri,r2) = Z^(r2,ri), that is U is symmetric under the exchange ri O r2. 

Let us consider this general expression in several different cases. Some of these cases look - only apparently 
singular: this behavior takes place as the triangle of the momenta ki degenerates into a line. 



• Equilateral case, 7*2 = r3 = 1. In this case we find 

Z^(l, 1) = Z^cquii = [6 ln(3/2) - 1] fc? w 1.43279 fc? . 
Since there is no apparent singular behavior for this result can also be confirmed numerically^. 

• Local case, r2 = 1, r3 — > (or r2 — > 0, r3 = 1). In this case, the limit exists and gives 

lim U (1, r-i) = Wiocai = -^fci ■ 

The result in the local case matches with that derived by Chen et al. [56| . 

• Singular line, = lime_>o(l ^ ^2 + e), with < r2 < 1. In this case we find 

„ r2{{r2 - 1)7'2[20 (r2 - 1) r2 + 9] + 6}(r2 - 1) - 6(r2 - 1)^ ln(l - r2) + Grf lnr2 



(B6) 



(B7) 



kl. 



24 (r2 - l)r2 

The limits r2 — ?> 0, 1 give again the local result, as expected. 
• Enfolded case, 7*2 = 1/2, = limc^o(l/2 + e). This is a particular case of the previous one. Then we find 

1 



(B8) 



limZi =Zionfoid = — (5 + 61n2)fc? w 0.38162 /c? 

e^O 24 



(B9) 



Mathematica, working in high precision, returns this same numerical value. The same value can also be found by applying the Method 
II mentioned above. 
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• Singular line, 7^3 = lime_^o(''2 ^ 1 + e), with r2 > 1. In this case we obtain 

(r2 - 1) {rsh (6r2 - 9r2 + 29) - 40] + 20} + Grf [(ra - l)^ ln(;:^) + Inrs 



lim 14 



kf. (BIO) 



e^o 24(r2-l)r2 
In the limit r2 — > 1^, this reproduces the local limit result, as expected. 
• Singular line, = lime_j.o(l + r2 — e). In this case we find 

r2{r2{r2[3r2(2r2 + 5) + 38] + 15} + 6} + 64(r2 + if + A + 6(r2 + 1)^^ ln(r2 + 1) 

limW= ^ ^ kf, (Bll) 

24r2(r2 + l) ' 

which recovers the local limit as r2 = e — > 0. 
Therefore, we have shown that the physical limits are all finite. 

The V integral 

The next step is to compute the integral ((8T|) . It is convenient to solve the integral by studying the limit 

ry 

V = lim lim / /v(xi)c?a;i + perm. , (B12) 

where /v is the integrand in Eq. (|8ip . As for the y — > 00 limit, we set the rapidly oscillating functions to vanish by 
regularizing e^''^^ as in the Method I of the integral U. After taking the limit e 0, we finally obtain 

V = iw ^ ^.J 3r2^ + Gr^j-I + rUrj - Qrg - 3) + r:,rl{6rl - 13r^ + 3) + 3rlri + 3 

3(^2 - ra - l)(r2 + rs - l)(r2 + r-s + 1) 

- Ibr^rl ~ Arjrl + 3r^(r3 - 1) - 12r|(r2 + 1) - r'2r3(19r| + Ibrj + ISr^ + 9) - 7r| + 3r| + 6r| + rj + Grg] 

6(7-2 - r3 + l)(r2 + r3 + 1) ^ ^ ^ ^ ^ ^ ^ ^ ^ 

- r|r2 + 6r|r2 + 5r3r2 + 19r2 + rl - - 6r| - 6r| + 5r3 + 7 - rlr2) 

kl In 2 



3(r2-r3-l)2(r2+r3-l) 
/c?ln(r2 +r3 + 1) 



2(^. ^rs- 1)2 ^^'^^ - 20ri - Gr^r^ + 22r2 + 4r2r2 - 20r2 + 4- Grj + 5) 

(4 + 2rzr\ - 2r\ - 3r^r^ - 2r3rf - 3?^ - 3rf r| + Sr^r^ - Srgr^ + 3r^ + 2r;^r2 - 2r|r2 



6 (r2 + 7-3 - 1) 2 

- 8r2r2 - 2r3r2 + 2r2 + - 2r;^ - 3r:^ + 3r2 + 2r3 - 1) 
fc3ln(l^) 

6 (7-2 - r3 - 1) ^ 

- 8r|r2 + 2r3?-2 + 2r2 - r| - 2r| + 3r| + 3r| - 2r3 - 1) 



{t% - 2r^r\ - 2r| - 3r2r| + 2r3ri - 3ri + 3r3r2 + 8r2r2 + Srgr^ + Sr^ + 2r4r2 + 2r|r2 



fc^ln 



''2 +''3 + 1 

2r3 



1 2 



{rl ~ 2t^t\ + 2r\ - irlrl ~ 2r^r\ - Zr\ + 3r|r2 - 8r2r2 + Srgrf - Srf + 2r4r2 - 2r37-2 



6 (7-2 -7-3 + 1) 

- 8r2r2 - 2r3r2 + 2r2 - r| + 2r| + 3r| - 3r2 - 2r3 + 1) , (B13) 
where V(r'i,r2) = V(r2,7'i). Let us now analyze this expression, on the lines/points of physical interest. 

• Equilateral case, 7-2 = 7-3 = 1. In this case we find 

V(l, 1) = Vequii = y [2 + ln(2/3)] kf « 11.959 kf . (B14) 

• Local case, r2 = 1, — ^ (or 7^2 — > 0, = 1). In this case, the limit exists and gives 

20 

lim V(l,r3) = Viocai = ^A:3, (B15) 
where this value of Viocai matches with the one derived by Chen et al. [56| . 
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Singular line, ^3 = lime_^o(l ^ ?'2 + e), with < r2 < 1. In this case we obtain 



{(r2-l)r2[ll(r2-l)r2-105]-86}r2(l-r2) + 6(r2-l)3ln(l-r2)-6r3lnr2 ,3 

hni V = — fci . (B16 

c^o 12 (r2 - 1) r2 

The limits r2 0, 1 again give the result (jB15| . 

• Enfolded case, r2 ~ 1/2, ~ lime_j.o(l/2 + e)- This is a particular case of the previous one. Then we find 

lim V = Vonfoid = f ^ - i In 2 ) kf « 4.5753 kf . (B17) 
«->o \ o4 z / 

• Singular line, — lime_j.o(?'2 — 1 + e), with r2 > 1. In this case we find 
(r2 - 1) [{r2[r2(86r2 - 105) + 94] + 22}r2 - 6(r2 - 1)^1 111(7^^) - ll] - Q4 In^s 



lim V 



kt . (B18) 



12 (r2 - 1) ra 

Taking the limit r2 1"*", we recover the result (jB15|) . 
• Singular line, = lime_>o(l + t'2 — e). In this case we obtain 

r2{r2{r2[r2{86r2 + 239) + 295] + 239} + 86) - 64{r2 + 1)-Mn( ;3_ + 1 ) _ 6(r2 + 1)^ ln(r2 + 1) 



hm V = — ^ '■ — z ^ — fc? , 

12r2 (r2 + 1) ' ' 

(B19) 

which again reproduces the value (jB15|) as r2 = e — ^ 0. 
Therefore, the integral V remains finite in the physical parameter space of r2 and r3. 
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